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ABSTRACT 

We examine the spectra and infrared colors of the cool methane-dominated atmospheres at T^ff < 1400 K ex- 
pected for young gas giant planets. We couple these spectral calculations to an updated version of the Marley 
et al. (2007) giant planet thermal evolution models that include formation by core accretion-gas capture. These 
relatively cool "young Jupiters" can be 1-6 magnitudes fainter than predicted by standard cooling tracks that in- 
clude a traditional initial condition, which may provide a diagnostic of formation. If correct, this would make true 
Jupiter-like planets much more difficult to detect at young ages than previously thought. Since Jupiter and Saturn 
are of distinctly super-solar composition, we examine emitted spectra for model planets at both solar metallicity 
and a metallicity of 5 times solar These metal-enhanced young Jupiters have lower pressure photospheres than 
field brown dwarfs of the same effective temperatures arising from both lower surface gravities and enhanced 
atmospheric opacity. We highlight several diagnostics for enhanced metallicity. A stronger CO absorption band 
at 4.5 /im for the warmest objects is predicted. At all temperatures, enhanced flux in K band is expected due to 
reduced collisional induced absorption by H2. This leads to correspondingly redder near infrared colors, which are 
redder than solar metallicity models with the same surface gravity by up to 0.7 inJ-K and 1. 5 in H-K. Molecular 
absorption band depths increase as well, most significantly for the coolest objects. We also qualitatively assess 
the changes to emitted spectra due to nonequilibrium chemistry. 

Subject headings: planets and satellites: formation; planetary systems; radiative transfer 



1. INTRODUCTION 

Astronomers around the world are making significant efforts 
to image planets in orbit around other stars CBeuzit et aLi2007t 
iNielsen et all l2007t iLafreniere etai] l2007t l^ai et all l2008l) . 
Work on suppressing the glare of potential parent stars has pro- 
ceeded to the point where contrast ratios of 10"^ can now typ- 
ically be achieved on the telescope, a nd 10"^ is on the horizon 
( [Macintosh et al.ll200d iDohlen et al.li2006,) . Since the contrast 
ratio for the Jupiter/Sun is 10" ^ (a contrast ratio which has now 
been achieved in a laboratorv.'Trau ger & Traubll2007l) . the ma- 
jority of this detection work focuses on young stars, as giant 
planets should be warmest, largest, and brightest when they are 
young, but will cool, contract, and f ade inexorably as they age 
joraboske et al. 1975; Bodenheimer 1976l;ISa"umon et al.lll996l; 
[Burrows et al.l ll997i) . Given the difficulty of these low contrast 
ratio observations, the interpretation of observed photometry 
and spectra takes on great importance. In practice, when faint 
planetary candidates are detected, evolution models, which aim 
to predict the structural and atmospheric properties with age, 
are needed to convert observed photometry or spectra into a 
probable planetary mass. 

The formation mechanisms of brown dwarfs and giant plan- 
ets are still not well understood in detail. While brown dwarfs 



likely form directly from molecular cl oud gas in something 
akin to the star form ation mechanism dLuhman et alj l2007al ; 
IWhitworfli et al.ll2007h . "true planets" form in a disk (lAU def- 
initions aside — for a discussion see Cha brier et al. 2007), per- 
haps predominantly via core accretion dLissauer & StevensonI 
l2007h . Recently, discussion has turned to how these distinct 
formation mechanisms, which may overlap at several Jupiter 
masses, may leave observational signatures in terms of an ob- 
ject's orbit, evolution, and atmosphere. A given parent star 
may well harbor both classes of low mass objects. This paper 
addresses the atmospheres of extrasolar giant planets (EGPs), 
while also describing the spectral properties of our recent work 
to couple core accretio n formation to subsequent planetary 
evolution (iMarlev et al.ll2007b ). In this introduction we will 
first review giant planet evolution models, then discuss our 
current understanding of the metal-enhanced atmospheres of 
Jupiter and Saturn. In §2 we describe our model atmosphere 
code, while in §3 we discuss the differences in atmospheric 
pressure-temperature (P-T) profiles, chemistry, and spectra be- 
tween models at solar metallicity and those at 5 x solar. §4 fo- 
cuses in particular on the near and mid infrared colors of metal- 
enhanced atmospheres while in §5 we discuss and tabulate the 
near- and mid-IR colors for our "hot start" and core-accretion 
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start evolution models. §6 addresses nonequilibrium chemistry 
while §7 contains additional discussion, caveats, and our con- 
clusions. 

1.1. The Early Evolution of Giant Planets 

Over the past decade only a small number of workers have 
attempted the difficult task of coupling non-gray radiative- 
convective atmosphere models to thermal evolution mod- 
els to enable an understanding of interior structure, atmo- 
spheric structure, atmospheric chemistry, and emitted spectra 
for giant planets and brown dwarfs (e.g.. Burrows et al. Il997t 



Chabrier et all l2000t iBaraffe et alJ 12003 
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200a) . It has perhaps only recently become appreciated by the 



wider community that these models do not include a mech- 
anism for the formation of the objects that they aim to un- 
derstand. The starting point for these models is an arbitrarily 
large and hot, non-rotating, adiabatic sphere. These model ob- 
jects are then allowed to cool and contract from this arbitrary 
state. The initial model is soon unimportant, as the cooling and 
contraction are initially very fast, since the Kelvin-Helmholtz 
time, fKH is inversely proportional to both luminosity and ra- 
dius. Although it is true that the models forget their initial 
conditions eventually, it is not immediately obvious how long 
this may take. In the past, a common thought was that af- 
ter "a few million years" the initial conditions are forgotten 
and that these standard "hot start" evolution models are reli- 
able. Although this type of m odel has been successfully ap- 
plied to Jupiter for decades (e.g., Grab oske et al.' 19751: iHubbardI 
ll977l:lGuillot et alJl995l:lFortnev & Hubbard 2003) their appli- 
cation to planets at very young ages could potentially be suspect 
JStevenson 1982). More recently Baraffe et al. (2002) have in- 
vestigated similar issues for brown dwarfs. 

In order to better understand the properties of gas giant plan- 
ets at young ages, in lMarlev et al.l (|2007b) we undertook an in- 
vestigation of the early evolution of giant planets, with initial 
properties given b y a state-of-the-art mo del of planet formation 
by core accretion dHubickvi et al.ll2005h . rather than the tradi- 
tional (but arbitrary) initial condition which we termed a "hot 
start." As shown in Figure[Tl the post-formation properties of 
these planets are surprising. The model planets started their 
lives smaller and colder than their hot start brethren. The core 
accretion start models were less luminous by factors of a few to 
100, and the initial conditions were not forgotten for timescales 
of tens of millions to one billion years. The reason for the sig- 
nificant difference l ies in t he treatment of gas ac cretion (see 
iMarlev et aLll2007bl) . In the lHubickvi et all j2005l) models the 
accreting gas arrives at nearly free fall velocity to a shock inter- 
face at the protoplanet. The shock radiation transfer is not fol- 
lowed directly, but a shock jump condition from Stab ler et al.l 
( Il980l) is employed; this accretion luminosity is entirely radi- 
ated away, leading to the prominent luminosity spike in these 
models during gas accretion. The gas that finally accretes onto 
the planet is therefore relatively cold, low entropy gas. 

It is therefore enticing to imagine that one could use the early 
luminosity, T^n, and surface gravity to determine the formation 
mechanism of a faint planetary-mass companion. This may be 
possible, but we caution th at the current generation of core ac- 
cretion formation models dHubickvi et al.l 120051 : lAlibert et al.l 
l2005al: llkoma et al.ll2000h are still only ID representations of 
a 3D process. A detailed look at radiation transfer in the for- 
mation shock, as well as incorporating multi-dimensional ac- 
cretion, should be undertaken before accurate luminosities of 



young planets can be confidently predicted. There is, however, 
another promising avenue for determining "planethood." While 
a brown dwarf-like companion and its parent star would be ex- 
pected to share common elemental abundances, the same may 
not be true of a companion that formed in a disk via core accre- 
tion. 

1 .2. The Atmospheres of "True " Giant Planets 

The Galileo and Cassini spacecrafts have unambiguously 
shown us that the atmospheres of Jupiter and Saturn, respec- 
tively, are enhanced in heavy elements relative to the Sun. The 
Galileo Entry Probe measured the abundances of oxygen, car- 
bon, nitrogen, sulfur, and various noble gases in the atmosphere 
of Jupiter Except for oxygen, a n enhancement of so- 
lar for each element was found jAtreva et al.ll2003l) . although 
the oxygen abundance det ermination may have been hi ndered 
by meteorological effects dShowman & IngersoIllll998l) . Sat- 
urn's atmosphere is enhanced i n carbon by factor of ~10, from 
an analysis of Cassini spectra jFlasar et al. 2005) , and in phos- 
phoru s by a factor of ^7, from ISO spectra tde Graauw et al.l 
1997: IVisscher & Feglevir2005b . If Jupiter, Saturn, and EGPs 
formed through a common mechanism, we can expect EGP at- 
mospheres to have high metallicities as well. 

How this atmospheric metallicity (and indeed the ratios of 
specific elements) may be set, as a function of planet mass, 
orbital distance, disk mass, disk metallicity, etc., is still open 
territory. In particular, the relative importance of processes 
that have enriched the atmospheres of Jupiter and Saturn is 
still unclear. These potentially include plan etesimal bombard- 
ment and accumulation during formation (Owenetal. 

iri999L 

iGautier et al . 2001a,b; Guillot & Gladman 2000; Alibert et alj 
2005^, erosion of the heavy element core dStevensoij 1 1 985b 



Guillot et aL 20041). di rect accretion of metal-rich disk gas 



(Guillot & Huesoll2006l). and chemical fractionation within the 
planet (iStevenson & Salpetej 1 9771: iLod ders 2(304). Clearly ob- 
servations of extrasolar giant planet will shed light on giant 
planet formation. Here we will undertake a first step at explor- 
ing how the spectra and colors of uniformly metal-enhanced 
atmospheres differ from strictly solar composition. 

Since^ diere is some evidence from both observa- 
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tions OChauvin et alJ I2005D and th eorv dBossI 120011: 
iKroupa & Bouviedl2003l: llda & Linll2004l) that the "planetary" 
and "stellar" formation modes may overlap at several Jupiter 
masses, it will be important to be able to decipher a formation 
mechanism based on observable properties. Some of this work 
has progressed on the orbital dynam ics of given companions 
or clas ses of companions, such as bvlRib as & Miralda-Escudi 
( l2007h . who found different eccentricity distributions for the 
radial velocity planets above and below 4 Mj. Additionally, 
transiting planet mass and radius determinations allow for cal- 
culations of bulk planet density, which shed light on th e internal 
abund ance of heavy e leme nts (G uillot et al. 2006; Fort ney et al.l 
'2007- Bui Tows et al.ll2007l: iBa-affe et al.ll200 8). A hallmark of 
Jupiter and Saturn is that they are enriched in heavy elements 
compare d to the Sun, which is known from planet structur e 
models dPodolak & CameronI [19741 ISaumon & Guillotll2004 . 
These heavy elements are partitioned betwe en a dense core and 
an en r ichment withi n the H/He en velope (Hubb ard & MarlevI 
[19891 iGuillotI [19991 ISaumon & Guillot 2004). All planets 
which form in disks around young stars are expected to be 
enriched in heavy elements due to these disks possessing both 
abundant gas and solids (.Pollack et al... 1996: .Ida & Lia,2004l) . 
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Another rewarding pathway for differentiating planets from 
low mass brown dwarfs is from direct characterization of their 
atmospheres. One could spectroscopically measure the heavy 
element abundances of the H/He envel ope directly, as has long 
been done in the solar system (e.g. iGautier & Owed 119891; 
lEncrenaj|2005h . If Nature is able to form two kinds of objects 
with an overlapping mass distribution, these distinct forma- 
tion modes may leave distinct observable atmospheric metallic- 
ities in the atmosph eres of these objects (Chabrier et al, 200 7; 
iMarlev et al.ll2007"ah . Therefore, two classes of planetary mass 
companions may be revealed by their emitted spectra. 

There are then two tasks to be completed. We first will 
investigate the differences in spectra and infrared colors be- 
tween models of solar composition (M/H=[0.0]) and those with 
a metallicity enhanced by a factor of five (M/H=[0.7]). This 
metallicity enhancement is similar to that of Jupiter and Saturn, 
but still far removed from the 30-40 x solar enhancement (at 
least in carbon) that has bee n measured for Uranus and Nep- 
tune (IGautier &Owenlll989 t). This investigation will be done 
at the low gravities and effective temperature most relevant for 
EGPs, (log g<4.3, 7;ft<1400 K), whereas field T-type brown 
dwarfs with similar T^ffS in general have surface gravities 10 
times larger. Since we restrict ourselves to Teff < 1400 K, we will 
necessarily be targeting cloud-free CH4-rich, rather than CO- 
rich, atmospheres. We focus on these relatively cool objects 
since a detection of CH4 together with an age estimate would 
significantly constrain a planetary candidate's mass, even tak- 
ing into account the uncertainties in evolution models. Later we 
will examine the evolution of infrared sp ectra and colors at 5 x 
solar metallicity, specifically for the Marl ev et alj j2007bl) evo- 
lution models for giant planets. As anticipated, the lower T^ff 
and radii for these models lead to dramatically fainter absolute 
magnitudes compared to hot start 1 x solar models. 

2. MODEL DESCRIPTION 

We employ a ID model atmosphere code that has been 
used for a variety of planetary and substellar objects. Re- 
cently it has been u s ed for brown dwarfs ( Marlev et all 2002 ; 



Saumonetal. 2006, 2007) and EGPs ( Fortnev et al.M2005b , 
2006- .Marlev et al. 2007a^; iFortnev & Marlev , 2007^ The ra- 



diative transfer method was developed by Toon et al. ( 19891) 
and has in the past been applied to Titan (McKav et al. ''1989'), 
Uranus ( Marley & McKay 1999), Gliese 229b (Marleyetal. 
119961) . and brown dwarfs in general (Burrows et al. jl997l) . 
We use the elemental abundance data of .Lodders 120031) and 
compute chemical equilibrium compositions at metallic ities 
of 1 X and 5 x so lar, following iFeglev & LoddersI d 19941) and 
iLodders & Fegle y (2002, 2006). The chemistry calculations 
include "rainout," where refractory species are depleted from 
the atm osphere due to their condensa tion into cloud decks 
(ILodde rs 1999; B urrows & Sharpll 19991) . The spectra of brown 
dwarfs and our solar system's giant planets can only be re- 
produced when chemistry calculations incorporate this process 
jFeglev & LoddersI 119941; iMarlev et alj 120021; iBurrows et al] 
120021). We use the correlated-k method for the tabulation of 
gaseous opacities (Goody etal. 1989); our extensive opacit y 
database is described in lFreedman. Marlev. & LoddersI (120081) . 
The model atmosphere code is used to compute radiative- 
convective equilibrium pressure-temperature (P-T) profiles and 



low resolution spectra. High resolution spectra are computed 
separately using a full line-by-line radiative transfer code, 
which utilizes the same chemistry and opacity database. Since 
we are modeling warm planets relatively far from their parent 
stars, here we ignore stellar insolation. 

Although we restrict our metal-enhanced planet atmospheres 
and comparative model spectra to Teff < 1400 K, we have com- 
puted 1 X metallicity models up to 2400 K for use in evolu- 
tionary calculations. This is the grid of model atmospheres 
that serves as the upper boundary condition for modeling 
the planets' thermal evolution. Since cloud opacit y is pre- 
dicted to affect the evolut ion of these planets ( Lun ine et al.l 
119891: IChabrier et alJl2000l) . we have elected to include it in 
the evolution atmosphere g rid. We use the cloud model of 
lAckerman & Marlev (1200 lb to describe the location, vertical 
distribution, and particle sizes of major cloud-forming species 
corundum (AI2O3), iron (Fe), and forsterite (Mg2Si04). We as- 
sume an /sed sedimentation efficiency p arameter of 2, which 
best matches obse rvations of L dwarfs dCushing et alj 120071 ; 
IMarlev et al.ll2008l) . We have generated a grid of cloudy model 
atmospheres from log ^=3.0 to 4.5 and reff=500 to 2400 K, sup- 
plemented with cloud-free models at reff<500 K. At low Teff, 
the refractory clouds reside very deep in t he atmosphere and 
negligibly effect the spectra and structure ( Saumon & Marlev! 
200 8). A more expans ive version of this grid is used in 
ISaumon & Marlev! (|2008) to compute the thermal evolution of 
brown dwarfs down to reff=50 K. Our previous pl anet evolu- 
tion calculations, presented by dMarlev et al.ll200'7bl) . neglected 
cloud opacity. 

The inclusion of cloud opacity into the atmosphere grids 
leads to some difference s in the cooling curve s from Figure[T] 
compared to those from iMarlev et al.l (l2007bh . Perhaps most 
notable is that the core-accretion start planets begin their evo- 
lution even colder here, because the cloud opacity closes off 
an atmospheric radiative zone at pressures of several to tens 
of bars, which has the effect of leaving a cooler photosphere 
for a given interior adiabat. In general, radiative zones can 
form because the relevant atmospheric opacities, particularly 
those of water and the pressure-induced opacity of H2, are 
strongly wavelength-dependent, and opacity windows can ap- 
pear around 1 to 2/^m at temperatures of 1000 to to 2,000K. 
The overlap of the Planck function with these windows can al- 
low local radiative transport of energy whereby the tempera- 
ture profile becomes less steep than an adiabat. As the local 
temperature continues to fall and the Planck function moves 
red-ward still, opacities increase again, closing the radiative 
window and the local temperature profile can again return to 
an adiabat. Such a two-layered convective structure was pre- 
dicted for Jupiter at pressures of severa l kilobar and tempera- 
tures near 2000 K bv lGuillot et all (11994) . The same effect was 
subsequently seen in models of brown dwarfs and warm giant 
planet atmospheres (Ma rlev et alj 119961; iBurrows et al.l 119971 ; 
Allard et al. 2001), but since such objects are warmer than 
Jupiter the corresponding radiative region is higher in the at- 
mosphere at lower pressur^ 

Details of th e calculation of the pl anetary evolution models 
can be found in IMarlev et al.l (12007 b'). All planets are assumed 
to be compos ed of pure H/He env elopes, with Y = 0.234 (the 
value used in lHubickvi et al.l [20051) overlaying a dense olivine 



^ Jupiter's deep radiative window at 1 /im originally noted by Guill ot et all 41994 ) was later found to likely be closed by highly pressure-broadened alkali opacity 
(Guillot et al. 2004; Freedman et al. 2008). The potential radiative region in the young giant planet models arises from the longer-wavelength near-infrared windows 
in water opacity (JHK) which are not as strongly affected by alkali opacity. At depth in Jupiter the water opacity windows are closed by the strong pressure-induced 
opacity of H2 at kilobar pressures. 
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core. The core masses range from 16 (for the 1 Mj planet) 
to 19 M© (for the 10 Mj planet). 

3. ATMOSPHERES AND SPECTRA AT SUPERSOLAR 
METALLICITY 

Before we generate spectra and colors for these evolutionary 
models, it is first worthwhile to examine the P-T conditions for 
these cool planetary atmospheres. In Figure|2]are shown a col- 
lection of cloud-free P-T profiles at a surface gravity of log g 
(cgs)=3.67, representative of a young 4 Mj planet. Solar com- 
position profiles at 1400, 1000, 600 K are shown in black, while 
profiles with 5 x solar metallicity are shown in red. Solid dots 
indicate where the local temperature equals the Teff, illustrat- 
ing the mean photospheric pressure. Chemical boundaries are 
shown in dashed curves and cloud condensation boundaries are 
dotted curves, using the same black/red color scheme. 

Quite clearly, high metallicity atmospheres are everywhere 
warmer at a given Teff. The higher gaseous opacity of these 
atmospheres leads to a photospheric pressure that is neces- 
sarily lower at a given T^s; one cannot see a s deeply into a 
higher opacity atmosphere (e.g., ISaumon eT al. 1994). This 
difference in pressure is a factor of '--^2-2. 5 for these 5x so- 
lar atmospheres. This lower photospheric pressure for higher 
metallicity atmospheres can manifest itself in lower contri- 
butions to the gaseous opacity for some species, specifically 
those that are pressure dependent. These include collisional 
induced absorption (CIA) by H2 molecules (iBorvsowl l2002l 
where the opacity per unit volume is dependent on the square 
of the pressure), and Na and K, which have strongly pres- 
sure broadened optical absorption lines (Bu rrows et al.l I2OO0I; 
AUard et al . 2003). As has been shown in detail elsewhere (e.g., 
Lodders & F eglev 2002), condensation and chemical bound- 
aries are non-trivial functions of metallicity. The higher the 
metallicity, the higher the temperature at which initial conden- 
sation will begin. However, the CO/CH4 and N2/NH3 equal- 
abundance curves move to cooler temperatures with increased 
metallicity dLodders & Feg lev 2002). 

The spectra of the profiles plotted in Figure |2] should show 
strong absorption due to CH4 and H2O, with CO also appearing 
in the warmer profiles and NH3 in the cooler profiles. Spectra 
for a collection of models are shown in Figure |3] The spectra 
for solar metallicity models are shown in black, while those for 
the 5x solar models are shown in colors. We can investigate 
these spectra for diagnostics of high metallicity. 

The most prominent difference in spectra between the "1 x" 
and "5x" models is a flux enhancement in K band for all 5x 
models. This brightened K band peak wa s previously shown for 
a "5 x" solar model at 7;ft=1200 K by ( Chabrier et alj|2007h in 
their discussion on differentiating planets from brown dwarfs. 
This flux enhancement is tied to the CIA opacity of H2, as it 
is substantially larger in K band compared to other near in- 
frared wavelengths. For lower pressures photospheric pressures 
this CIA absorption is weaker, letting more flux escape in K 
band, relative to other wavelengths. Also apparent is the in- 
creased CO absorption at 4.5 /im. However, the greater CO 
absorption at 2.3 /im is swamped by the higher /T-band flux. At 
higher metallicity the abundances of metal-metal species such 
as CO are increased to a larger degree than those of hydride 
species. This leads to stronger CO absorption, relative to ab- 
sorption bands from other molecules. This 4.5 /im CO absorp- 
tion band will likely be a valuable diagnostic for the warmest 
young Jupiters (.Chabrier et al...2007.) . but its importance nec- 



essarily wanes as the CO abundance drops dramatically below 
Teff-^^IOOO K for these objects. Absorption band depths due to 
H2O and CH4 are modestly deeper at higher metallicity as well. 
At the coolest TeffS the bands are so deep that it may well be 
difficult to see any emitted flux, making band depths a difficult 
diagnostic in practice. 

The contrast ratios that will be needed to directly image can- 
didate planets are show in FigureH) The flux density at 10 pc 
is plotted for three models of a 4 Mj planet at an age of ~10 
Myr. Two models utilize the hot start initial condition (and al- 
low a comparison of the effects of metallicity only) while one 
uses the core accretion initial condition and also has a 5 x so- 
lar atmosphere. Solid curves show the necessary contrast ra- 
tios around a Sun-like star, while dashed curves are contrast ra- 
tios around an M dwarf. Strikingly, while the hot start models 
would be easily detectable at 10"^ contrast in the near infrared, 
the core-accretion start model would be undetectable. As men- 
tioned above, the H2O and CH4 bands are exceedingly dark, 
especially for the cooler high metallicity planets, such that con- 
trast ratios of 10"' may be necessary to see flux from inside 
these bands. 

Figure|5]allows one to look forward in time to examine these 
planets at age of 80 Myr By this age the two different cool- 
ing tracks have nearly run together, meaning the core-accretion 
model has nearly forgotten its initial condition. Although the 
contrast ratios for all bands in the near infrared are smaller for 
all the models, compared to 10 Myr, one now has the important 
and distinct advantage that the uncertainties due to the forma- 
tion history are much smaller. This will enable a more realistic 
mass estimate for the planet once one obtains photometric or 
spectral data. 

It is also well known that surface gravity changes pho- 
tospheric pressures in substellar and planeta ry objects (e.g., 
'Burr ows et all 1 19971 120021: iKirkpatricklllOOl . Since a given 
optical depth is proportional to column density (as long as CIA 
opacity does not dominate), this optical depth is reached at a 
higher pressure in high gravity objects, meaning they will have 
high pressure photospheres. This effect is illustrated in Figure |6] 
for solar metallicity atmospheres at 1200 and 600 K, at gravi- 
ties of log g='i.5, 4.5, and 5.5. Whereas old massive brown 
dwarfs commonly have gravities from 5.0 to 5.5, young Jupiters 
likely exhibit gravities from 3.5 to 4.0, meaning they will have 
warmer atmospheres at a given T^s- Admittedly, Figure |6] plots 
relative extremes in gravity, but the large differences in photo- 
spheric press ures are significant. A s has been shown by other 
authors (e.g.. iBurrows et aT]|2003h lower gravity atmospheres 
must cool to lower Teff values before crossing a given chemi- 
cal or condensation boundary. While the log g=5.5 model has 
clearly passes the water cloud condensation curve at 600 K, the 
log g=3.5 model does not. 

It is well understood that lower gravity and higher metal- 
licity both lead to lower pressure photospheres. Since young 
Jupiter atmospheres couple both higher metallicity and lower 
surface gravity it may well be that the visible effects of, for in- 
stance, the L-to-T transition, disequilibrium chemistry due to 
vertical mixing, and water cloud condensation may differ in de- 
tails compared to high gravity brown dwarfs. In §6 we quali- 
tatively address the effects of vertical mixing on carbon chem- 
istry. Differences in particles sizes, opacity, and vertical extent 
of clouds are possible, even likely. Simple cloud models predict 
cloud particles sizes t han may be ^10 times larger in these low 
gravity environments jAckerman & Marlevl2001l;ICooper et al.l 
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I2003h which may lessen the opacity of cloud decks. While it is 
still not clear what causes the rapid drop in the opacity of the sil- 
icate and iron clouds at the L-to-T transition in brown dwarfs, 
it should be kept in mind that this transition could be differ- 
ent in character at surface gravities up to 100 times smaller. 
For instance, there is already some evidence that the L-to-T 
transition may occur at lower Teg in low gravity b rown dwarfs 
jMetchev & Hillenbra nd 2006: Luh man et al.l '2007b), which is 
reasonable since lower gravity atmospher es trace the P-T space 
of higher gravity objects at a lower T^ff dBurrows et alj|2002t 
iKnapp et alJl2004 . 

4. "METALLICITY COLOR": QUANTIFYING METALLICITY 
EFFECTS IN INFRARED BANDS 

In Figure[3] we showed spectra of a subset of our 1 x and 
5x solar models. We can further investigate the differences 
between these models in the commonly used near and mid IR 
bands. Independent of an evolution model (since these models 
are uncertain at young ages) we can compare the differences in 
magnitudes at various values of Teff and log g, assuming no dif- 
ference in radii between 1 x and 5 x solar modelfl Here we 
will term the difference in magnitude between the 5 x and 1 x 
solar models the "metallicity color," MC=mB5-mBi where the 
letter subscript B is the band and the numerical subscript is the 
metallicity. A negative MC means that flux is enhanced in the 
higher metallicity model. 

Figure |2] shows the MC as function of Teff at three values 
of log g that span the gravity values for the models shown in 
Figurel?] Although this is an unusual way to plot colors, it al- 
lows for a clear determination of the best diagnostics for high 
metallicity. From 1400 K down to 800 K, J-K colors are --0.7 
redder for 5 x models than for 1 x models. Below 700 K higher 
metallicity depresses H band flux and even more greatly en- 
hances K band flux, leading Xo H-K colors becoming even 
redder than the J-K colors, with H-K high metallicity mod- 
els becoming redder by 2.0-2.5 than solar composit ion models, 
although water cl oud condensation below 500 K ( Marl ev et al.l 
l2002l:lBurrowset al. 2003) could modify this steep falloff. The 
L' band is interesting in that above 800-900 K flux is enhanced 
in L' at increased metallicity (see Figure[3]) but the flux is de- 
pressed in L' for the high metallicity models below reff'^800 
K. The remaining bands, Y, Z, and M' are marginally dimmer 
in the higher metallicity models, with a flux depression that is 
monotonically weaker for bands Y and Z as gravity increases. 

5. MAGNITUDES AND COLORS FOR YOUNG JUPITER 
EVOLUTION MODELS 

W e have coupled our spectral models to our updated 
iMa rlev et al. (2007b) evolution models show in Figure[T] This 
aUows for a determination of absolute magnitudes for these 
model planets as a function of age. In Table [T] we provide a ta- 
ble of Teff, radius, as well as absolute magnitudes in the standard 
red-optical and near infrared filters for the hot start 1 x models, 
which include cloud opacity. In Table|2]we provide a similar 
table for the core-accretion start models, assuming cloud-free 
5 X solar metallicity. This coupling for the core-accretion start 
models is not strictly self-consistent as the original evolution 
models used 1 x atmospheres, but this difference is a small one 
compared to the uncertainty in the initial condition. Since initial 
TeffS are all below ^800 K, the core-accretion start models span 

"* However, a planet with a 5 X solar abundance of metals mixed throughout the entire H/He envelope would be modestly smaller. The effect of a "solar" amount of 
envelope metals on the structure of brown dwarfs is commonly ignored tChabrier & Baraffe..2000, ; ,Saumon & Marlev.,2QQ8i) . 



a much smaller range of reff than our own hot start models and 
those of other authors (e.g.. iBurrows et al]|l997l : [Baraffe et al.l 
I200I . Here we also limit our calculations to T^a > 400 K, since 
we also ignore opacity due to water clou ds, which condense 
around reff=400-500 K for these objects jMarlev et all 120021: 
iBurrows et ani2003l) . We note that lBurrows et alj (120031) have 
shown that water cloud opacity has a relatively modest effect 
on the spectra of very cool objects. We will also address this 
issue in a forthcoming work on the spectra of ultracool dwarfs 

and p lanets below 7 K. 

In IMarlev et al] (12007b') and iFortnev et all (l2005ah we 
showed that evolution models that incorporate a core accre- 
tion initial condition are significantl y less luminous tha t our 
own "hot start" mod els and those of iBaraffe et al.l (l2003l) and 
IBurrows et alj (Il997h . For example, the luminosity difference 
for a 4 Mj planet is initially a factor of ~ 100 at 1 Myr (See 
Figure[T]) Given this factor of 100, one would expect a differ- 
ence in absolute magnitude of five. This is indeed born out, as 
is shown in Figure [8] which compares absolute magnitudes for 
a 4 Mj planet from our hot start 1 x and core-accretion start 5 x 
calculations. As indicated in Figure[T] even by ^100 Myr the 
hot start and core accretion cooling tracks have not yet merged. 
This behavior is apparent as well in Figure[8] which show ^0.5- 
1.5 magnitude differences at 100 Myr 

6. EFFECTS OF NON-EQUILIBRIUM CHEMISTRY 

Another important detection issue is in which bands searches 
for young Jupiters would be most efficiently carried out. From 
1400 K down to 500 K, the peak in planeta ry flux gradu- 
ally shifts from L' to M' as CO is lost to CH4 ( Burrows et al.l 
1997). The 4 Mj core-accretion planet, which is always rel- 
atively cool, is two magnitudes brighter in M' band than L' . 
However, the M'-band flux will likely be d epressed somewhat 
due to non-equilibriurn CO/CH4 chemistry dFegley & LoddersI 
1996: Nolletal. 1997: Saum on et al.ll2003L l2006h. Rece"ntiv 
Hubeny & Burrows (2007) have investigated these affects for 
solar composition models across a wide range of Teff, but at 
higher gravity than we consider here, and find a ^40% flux 
decrement in M' band due to non-equilibrium chemistry. In- 
terestingly, they find that this flux decrement increases with de- 
creasing gravity (for a given K,,, the eddy diffusion coefficient), 
meaning that non-equilibrium chemistry will likely remain im- 
portant for young Jupiters. Indeed, the effects of vertical mixing 
on CH4/CO mixing ratio s were first descri b ed in Jupiter itself 
dPrinn& Barshavi 119771: lYung et alj [19881 iFeglev & LoddersI 
11994 . Although planets may well remain brighter in M' than 
L', the lower thermal background in L' for ground-based obser- 
vatories will make both bands attractive for EGP searches. 

While chemical reaction time constants can be measured (al- 
though there is uncertainly as to the actual reaction pathway 
for carbon chemistry), the mixing time scales are much more 
uncertain. In the convective atmosphere the mixing time scale 
can be computed from mixing length theory. However, in the 
radiative region, there is no a priori theory to characterize this 
mixing, parametrized by X.,. Previous studies in brown dwarfs 
have varied this coefficient over 6-8 orders of magnitude. Para- 
metric studies of the effects of non-equilibrium chemistry in 
brown d warfs, as as function of r^ff, gravity, and K^^ can b e 
found in lSaumonetal. (2003) and lHubeny & Burrowsl (l2007l) . 
In some instances K^ can be constrained by observations of 
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brown dwarfs (e. g. I Saumon et alj|2006h . There is also a long 
history of modeling non-equiUbrium chemistry in Jupiter, and 
recent estimates of K^^ in Jupiter's radiative atmosphere range 
from 10^ to 10"* (Bezard et al. 2002; Moses et al. 2005). There 
has not yet been a study of non-equilibrium chemistry at the low 
gravities relevant for young gas giants. When abundant data for 
these planets becomes available, detailed studies will of course 
be necessary. For now we will treat the effects of nonequilib- 
rium chemistry as an uncertainty in the model spectra, which 
we can briefly gage. 

Us ing the methods described in I Saumon et all (120031 l2006l 
l2007h we have computed models with K^^ = 10"* cm^ s"\ for 1 x 
and 5x solar metallicity, which we show in Figure |9] These 
near-infrared spectra clearly show the effects of enhanced CO 
and depleted CH4. Absorption by CO in M' band remains 
strong down to 500 K, leading to a large flux decrement, es- 
pecially in the 5x solar model. In the H, K, and L' bands, 
the reduced CH4 abundances leads to greater flux (sometimes 
dramatically) escaping through these bands. The behavior as 
a function of gravity, metallicity, and is complex due to a 
number of factors related to the regions of P-T space that ob- 
servations are sensitive to. For instance, lower gravity objects 
have lower pressure photospheres at a given T^ff, while higher 
atmospheric metallicity pushes the CO/CH4 transition to higher 
pressure, while higher values of lead to gas being mixed up 
from higher pressures and higher temperatures. 

We caution that these nonequilibrium spectra should only be 
viewed as illustrative. The non-equilibrium spectra are com- 
puted with the same P-T structure as the equilibrium models, 
even though the gas composition and the opacity have changed. 
The resulting non-equilibrium spectra have a larger integrated 
flux than the corresponding equilibrium spectra. In the fu- 
ture we will compute P-T profiles that are consistent with the 
nonequilibrium chemical abundances, alleviating this problem 
(e. g. ,Hubenv & Burrows. .2007). 

7. DISCUSSION & CONCLUSIONS 

We have shown how metal enhanced atmospheres differ from 
their solar composition counterparts, in atmospheric structure, 
chemistry, spectra, and colors. We have app l ied thes e re- 
sults to an updated version of the iMarlev et alj ( l2007bl) evo- 
lution models, which give cooling tracks for EGPs that are 
initially significantly colder than traditional models. We urge 
caution in the application of the computed absolute magni- 
tudes for these mod els provided in Table |2l Recall that the 
iMa rlev et al) (l20 07b|) mo dels incorporate the formation models 
ofTHubickyj et al.i (l2005h . which employ a treatment of accre- 
tion that is surely much simpler than what occurs in Nature. The 
potential agreement or disagreement between observations and 
the model cooling tracks and magnitudes should not be taken 
as evidence for or against the viability of the c o re acc retion for- 
mation scenario. Indeed, the iHubickvi et alj ( l2005h prescrip- 
tion is just one of several models of core accretion, which all 
currently include simplifications of the gas and solid accretion. 
While we can claim with some confidence that young Jupiters 
are fainter than those predicted from an arbitrarily hot start, 
how much fainter depend s sensitively on the details of accre- 
tion jMarlev et alj 1200 75"). Given the difficulty of predicting 
properties of EGPs at young ages, observations of these young 
objects will be of central importance. 

The next generation of direct imaging platforms will 
be the Gemini Planet Imager (GPI), likely at Gem- 



ini South ( [Macintosh et al.l l2006l) . and Spectro-Polarimetric 
High-contrast Exop lanet REsearch (SPHERE), at the VLT 
dDoWen et al.' 2006). The GPI will target the YJHK bands, 
while SPHERE will focus on YJH. Both instruments will in 
particular emphasize H band. In Figure[TO] we show a four 
panel plot that illustrates changes to //-band fluxes due to T^ff, 
gravity, metallicity, and nonequilibrium chemistry. All are ref- 
erenced from a standard case model with reff=700 K. Since GPI 
and SPHERE will likely employ custom narrow-band filters, 
we have overlain in gray possible narrow filters, kindly pro- 
vided by J. Graham. It is clear that the more limited the ob- 
servations are in wavelength, the more difficult planetary char- 
acterization will be, as lower surface gravity and higher metal- 
licity generally effect spectra in similar ways. We must also 
caution that although our understandin g of the gaseous opac- 
ity in these atmospheres is improving (Sharp & Burrowsl l2007l : 
Freedman et al. 2008), calculations of the contribution due to 
CH4 are quite uncertain. The absorption cross-section of CH4 
is difficult to model under the relevant P-T conditions found 
here, which manifests itself in mismatches of our models to 
brown dwarf spectra, especially around //-band dSaumon et al.l 
2006, 20071). 

Caveats now aside, we can readily summarize our findings 
for the low-gravity metal-enhanced young Jupiter atmospheres 
into the following four points: 

1. Young Jupiter atmospheres will have lower pressure 
photospheres than old field brown dwarfs due to their 
lower surface gravity (which has long been understood) 
and higher atmospheric opacity (if the planets have high 
atmospheric metallicities, like Jupiter and Saturn). 

2. Higher metallicity atmospheres, while generally hav- 
ing more opacity at all wavelengths, have relatively less 
opacity in /T-band relative to other bands, due to weak- 
ened CIA H2 opacity. This leads to a /T-band flux en- 
hancement of ^0.5 to 1.0 magnitudes between TeffS of 
500 to 1400 K. 

3. A spectral signature of high metallicity at Teff > 1000 K 
is a markedly deeper CO absorption band at 4.5 /im. 

4. A photometric feature of high metallicity at reff< 1400 
K is redder J-K and H-K colors, which may be redder 
by -0.7-1.5. 

We note that points 1-3 in particular echo the findings of 
IChabrier et al.l (12007 ). who had previously analyzed the spec- 
trum of a representative metal-enhanced planet model within 
the parameter range we examine here. The agreement is en- 
couraging. 

In closing, we note that the current best example of how 
well we may eventually be able to constrain the properties of 
a young EGP comes from the cool T7.5 dwarf Gliese 570D. 
This brown dwarf is a wide companion to the well-studied K4 
V stai- Gl 570A, and to a pair of M dwarfs, Gl 570BC. As dis- 
cussed by Saumon et al.l (2006), the distance, metallicity, and 
age of the system effectively constrain the physical parameters 
of the T dwarf Gl 570D. The spectrum is extremely well sam- 
pled from visible wavelengths, across the nea r infrared, to the 
mid infrared with Spitzer IRS. ISaumon et aL |(|2006) constram 
reff=800-820K, log ^=5.09-5.23, log (L/L(:i)=5.525-5.551, and 
mass=38-47 Mj. ISaumon et al.l (|2007|) additionally investigated 
two late field T dwarfs with similar spectral coverage, but with- 
out parent stars to constrain metallicities. Uncertainties in Tgff 
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increase ~50-100 K, while mass estimates widen to M «30- 
60 Mj. Since young Jupiters will have unknown metallicities 
and early on these objects will sample gravity and metallicity 
ranges that only marginally overlap the more well understood 
brown dwarfs, it will be challenging to constrain planetary pa- 
rameters with limited photometric and spectral data. 
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Table 1 

1 X Solar Absolute Magnitudes for "Hot Start" Evolution Models 



Mass (Mj) 


t (yrs) 


TLrr (K) 


R(R,) 


My 


Mz 


Mj 


M„ 


Mk 


M^, 


M„, 


LO 


l.OOOe+06 


900.3 


1.57 


16.20 


16.97 


15.05 


14.44 


14.01 


12.52 


11.22 




2.154e+06 


747.5 


1.53 


17.28 


18.13 


16.09 


15.67 


15.14 


13.49 


11.58 




4.642e+06 


644.2 


1.43 


17.82 


18.80 


16.75 


16.87 


16.37 


14.31 


12.04 




l.OOOe+07 


554.6 


1.35 


18.22 


19.35 


17.41 


18.27 


17.87 


15.19 


12.63 


2.0 


l.OOOe+06 


1266.7 


1.72 


14.36 


15.00 


13.18 


12.41 


1 1.72 


10.43 


10.45 




2.154e+06 


1048.8 


1.57 


15.36 


16.05 


14.23 


13.50 


13.12 


1 1.69 


10.99 




4.642e+06 


855.3 


1.47 


16.64 


17.42 


15.46 


14.84 


14.40 


12.87 


1 1 .40 




1 .0()0e+07 


710.3 


1.39 


17.53 


18.41 


16.36 


16.05 


15.61 


13.80 


11.83 




2.154e+07 


605.8 


1.33 


17.90 


18.96 


16.94 


17.36 


17.04 


14.65 


12.35 




4.642e+07 


512.9 


1.27 


18.42 


19.59 


17.71 


18.94 


18.85 


15.63 


13.02 


4.0 


1 .OOOe+06 


1657.0 


1.92 


12.66 


13.29 


11.48 


10.91 


10.19 


9.18 


9.43 




2.154e+06 


1432.0 


1.71 


13.73 


14.35 


12.53 


11.84 


11.09 


9.92 


10.15 




4.642e+06 


1207.6 


1.56 


14.79 


15.43 


13.62 


12.80 


12.28 


10.97 


10.74 




l.OOOe+07 


989.1 


1.45 


15.88 


16.58 


14.71 


13.95 


13.60 


12.14 


11.16 




2.154e+07 


805.3 


1.38 


17.00 


17.80 


15.80 


15.20 


14.84 


13.16 


11.57 




4.642e+07 


673.8 


1.32 


17.73 


18.66 


16.59 


16.43 


16.18 


14.04 


12.07 




l.OOOe+08 


574.6 


1.26 


18.04 


19.13 


17.09 


17.75 


17.80 


14.93 


12.66 


6.0 


l.OOOe+06 


1984.0 


2.10 


11.58 


12.19 


10.50 


10.01 


9.40 


8.58 


8.75 




2.154e+06 


1720.9 


1.82 


12.63 


13.25 


11.46 


10.89 


10.22 


9.25 


9.45 




4.642e+06 


1464.2 


1.63 


13.70 


14.31 


12.50 


11.81 


11.10 


9.97 


10.19 




l.OOOe+07 


1219.6 


1.50 


14.83 


15.45 


13.65 


12.82 


12.32 


11.04 


10.79 




2.154e+07 


988.4 


1.40 


16.02 


16.72 


14.81 


14.03 


13.68 


12.21 


11.22 




4.642e+07 


803.5 


1.33 


17.06 


17.86 


15.80 


15.25 


14.91 


13.19 


11.63 




l.OOOe+08 


674.9 


1.28 


17.81 


18.73 


16.58 


16.45 


16.25 


14.04 


12.13 




2.154e+08 


574.9 


1.22 


18.12 


19.19 


17.07 


17.75 


17.89 


14.94 


12.73 


8.0 


l.OOOe+06 


2184.9 


2.28 


10.91 


11.48 


9.92 


9.44 


8.89 


8.17 


8.30 




2.154e+06 


1987.6 


1.93 


11.74 


12.34 


10.66 


10.18 


9.60 


8.78 


8.97 




4.642e+06 


1656.0 


1.69 


12.91 


13.52 


11.72 


11.14 


10.47 


9.48 


9.70 




l.OOOe+07 


1400.7 


1.54 


14.08 


14.68 


12.87 


12.13 


11.44 


10.29 


10.43 




2.154e+07 


1163.4 


1.42 


15.15 


15.78 


13.97 


13.14 


12.71 


11.38 


10.95 




4.642e+07 


940.0 


1.34 


16.39 


17.10 


15.13 


14.39 


14.03 


12.51 


11.37 




l.OOOe+08 


769.1 


1.28 


17.35 


18.16 


16.04 


15.59 


15.29 


13.43 


11.81 




2.154e+08 


659.2 


1.22 


17.93 


18.86 


16.66 


16.69 


16.59 


14.21 


12.29 




4.642e+08 


543.2 


1.17 


18.41 


19.49 


17.35 


18.23 


18.59 


15.28 


13.02 


10.0 


1.0()0e+06 


2315.7 


2.44 


10.47 


11.00 


9.54 


9.07 


8.55 


7.89 


7.98 




2.154e+06 


2168.4 


2.04 


11.18 


11.75 


10.18 


9.71 


9.17 


8.45 


8.61 




4.642e+06 


1873.0 


1.75 


12.26 


12.87 


11.14 


10.63 


10.03 


9.15 


9.34 




l.OOOe+07 


1553.8 


1.57 


13.34 


13.94 


12.15 


11.52 


10.85 


9.80 


10.06 




2.154e+07 


1307.0 


1.44 


14.58 


15.19 


13.38 


12.58 


11.98 


10.77 


10.72 




4.642e+07 


1072.7 


1.35 


15.68 


16.33 


14.47 


13.66 


13.30 


11.89 


11.17 




1.0()0e+08 


871.5 


1.28 


16.83 


17.58 


15.54 


14.88 


14.55 


12.89 


11.57 




2.154e+08 


735.9 


1.22 


17.60 


18.45 


16.29 


15.95 


15.76 


13.68 


12.02 




4.642e+08 


614.5 


1.16 


18.07 


19.08 


16.88 


17.24 


17.43 


14.60 


12.61 




l.OOOe+09 


498.2 


1.12 


18.88 


19.97 


17.84 


18.82 


19.51 


15.67 


13.36 



Note. — Atmospheric metaUicity is 1 x solar. MKO filter set used. Time steps in years are equally spaced in log t. 
Only models with Tcf!> 500 are tabulated. 
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Table 2 

5 X Solar Absolute Magnitudes for Core-Accretion Start Evolution Models 



Mass (Aij) 


t (yrs) 




K (Kj) 


Wij 


My 




My 


M/y 


Ma' 


Mi/ 


M„/ 


Mc//4 


1.0 


l.OOOe+06 


672.5 


1.46 


21.70 


16.75 


18.03 


16.16 


17.07 


16.15 


14.63 


12.19 


16.23 




2.154e+06 


641.6 


1.43 


22.03 


17.00 


18.31 


16.44 


17.50 


16.52 


14.89 


12.33 


16.59 




4.642e+06 


592.9 


1.38 


22.58 


17.43 


18.76 


16.91 


18.22 


17.12 


15.32 


12.58 


17.21 




l.OOOe+07 


528.6 


1.32 


23.51 


18.13 


19.48 


17.72 


19.43 


18.07 


15.94 


12.98 


18.30 




2.154e+07 


455.8 


1.26 


24.76 


19.21 


20.56 


18.95 


21.05 


19.40 


16.81 


13.55 


19.78 


2.0 


l.OOOe+06 


652.2 


1.36 


22.00 


16.97 


18.27 


16.36 


17.40 


16.47 


14.84 


12.39 


16.46 




2.154e+06 


641.1 


1.35 


22.11 


17.06 


18.36 


16.45 


17.55 


16.60 


14.92 


12.43 


16.58 




4.642e+06 


625.7 


1.34 


22.26 


17.18 


18.49 


16.59 


17.76 


16.77 


15.04 


12.50 


16.76 




l.OOOe+07 


599.6 


1.32 


22.52 


17.39 


18.71 


16.82 


18.11 


17.06 


15.25 


12.62 


17.05 




2.154e+07 


550.3 


1.29 


23.16 


17.90 


19.23 


17.40 


18.98 


17.76 


15.71 


12.92 


17.84 




4.642e+07 


483.5 


1.25 


24.09 


18.69 


20.02 


18.30 


20.26 


18.81 


16.41 


13.36 


18.99 




l.OOOe+08 


409.0 


1.21 


25.38 


19.89 


21.19 


19.65 


22.06 


20.37 


17.40 


13.99 


20.63 


4.0 


l.OOOe+06 


585.3 


1.26 


22.71 


17.57 


18.87 


16.97 


18.32 


17.29 


15.36 


12.78 


17.19 




2.154e+06 


584.2 


1.26 


22.73 


17.58 


18.88 


16.98 


18.34 


17.30 


15.37 


12.78 


17.20 




4.642e+06 


580.6 


1.26 


22.77 


17.62 


18.92 


17.02 


18.40 


17.35 


15.41 


12.81 


17.26 




l.OOOe+07 


573.1 


1.26 


22.86 


17.69 


19.00 


17.11 


18.52 


17.46 


15.48 


12.85 


17.37 




2.154e+07 


558.8 


1.25 


23.03 


17.84 


19.14 


17.27 


18.76 


17.66 


15.61 


12.94 


17.58 




4.642e+07 


536.7 


1.24 


23.29 


18.06 


19.37 


17.52 


19.13 


17.97 


15.82 


13.07 


17.91 




l.OOOe+08 


492.8 


1.22 


23.83 


18.54 


19.84 


18.05 


19.90 


18.63 


16.26 


13.35 


18.60 




2.154e+08 


428.1 


1.19 


24.84 


19.52 


20.80 


19.14 


21.43 


20.02 


17.12 


13.89 


19.98 


6.0 


1 .OOOe+06 


563.5 


1.21 


22.98 


17.83 


19.12 


17.22 


18.64 


17.63 


15.58 


12.97 


17.44 




2.154e+06 


562.9 


1.21 


22.99 


17.84 


19.13 


17.22 


18.65 


17.64 


15.58 


12.98 


17.45 




4.642e+06 


561.2 


1.21 


23.01 


17.85 


19.15 


17.24 


18.68 


17.66 


15.60 


12.99 


17.47 




l.OOOe+07 


557.7 


1.21 


23.05 


17.89 


19.18 


17.28 


18.73 


17.71 


15.63 


13.01 


17.52 




2.154e+07 


550.9 


1.21 


23.12 


17.96 


19.25 


17.36 


18.84 


17.81 


15.69 


13.05 


17.62 




4.642e+07 


539.2 


1.20 


23.25 


18.07 


19.36 


17.48 


19.03 


17.97 


15.80 


13.12 


17.78 




l.OOOe+08 


517.2 


1.19 


23.50 


18.29 


19.58 


17.73 


19.38 


18.28 


16.01 


13.25 


18.09 




2.154e+08 


477.9 


1.18 


24.02 


18.79 


20.07 


18.27 


20.17 


19.01 


16.47 


13.54 


18.80 




4.642e+08 


423.3 


1.15 


24.83 


19.59 


20.85 


19.16 


21.45 


20.21 


17.21 


14.00 


19.95 


8.0 


l.OOOe+06 


556.3 


1.17 


23.09 


17.95 


19.23 


17.32 


18.74 


17.78 


15.67 


13.08 


17.50 




2.154e+06 


556.0 


1.17 


23.09 


17.95 


19.24 


17.32 


18.74 


17.78 


15.67 


13.08 


17.50 




4.642e+06 


555.2 


1.17 


23.10 


17.96 


19.24 


17.33 


18.75 


17.80 


15.68 


13.09 


17.51 




l.OOOe+07 


553.2 


1.17 


23.12 


17.98 


19.26 


17.35 


18.79 


17.82 


15.70 


13.10 


17.54 




2.154e+07 


549.4 


1.17 


23.16 


18.02 


19.30 


17.39 


18.85 


17.88 


15.73 


13.12 


17.59 




4.642e+07 


541.2 


1.17 


23.25 


18.10 


19.38 


17.48 


18.97 


17.99 


15.81 


13.17 


17.71 




l.OOOe+08 


524.4 


1.16 


23.43 


18.26 


19.55 


17.66 


19.24 


18.23 


15.97 


13.27 


17.94 




2.154e+08 


496.5 


1.15 


23.74 


18.55 


19.83 


17.98 


19.70 


18.65 


16.24 


13.45 


18.35 




4.642e+08 


462.1 


1.13 


24.23 


19.05 


20.31 


18.52 


20.50 


19.42 


16.71 


13.74 


19.07 


10.0 


1 .OOOe+06 


561.0 


1.14 


23.07 


17.96 


19.23 


17.29 


18.66 


17.76 


15.65 


13.12 


17.42 




2.154e+06 


560.8 


1.14 


23.07 


17.96 


19.23 


17.30 


18.67 


17.77 


15.66 


13.12 


17.42 




4.642e+06 


560.4 


1.14 


23.08 


17.96 


19.24 


17.30 


18.67 


17.77 


15.66 


13.12 


17.42 




l.OOOe+07 


559.0 


1.14 


23.09 


17.98 


19.25 


17.31 


18.69 


17.79 


15.67 


13.13 


17.44 




2.154e+07 


556.1 


1.14 


23.12 


18.01 


19.28 


17.34 


18.74 


17.83 


15.70 


13.15 


17.48 




4.642e+07 


550.0 


1.14 


23.18 


18.06 


19.34 


17.41 


18.83 


17.92 


15.76 


13.18 


17.56 




l.OOOe+08 


537.8 


1.13 


23.31 


18.18 


19.46 


17.54 


19.02 


18.09 


15.87 


13.26 


17.73 




2.154e+08 


518.2 


1.12 


23.51 


18.37 


19.65 


17.75 


19.32 


18.37 


16.05 


13.38 


17.99 




4.642e+08 


489.7 


1.11 


23.84 


18.69 


19.96 


18.10 


19.83 


18.86 


16.36 


13.57 


18.45 




l.OOOe+09 


442.1 


1.09 


24.49 


19.36 


20.61 


18.83 


20.93 


19.94 


17.01 


13.98 


19.43 



Note. — Atmos pheri c metallicity is5x solar. MKO filter set used. Filter "CH4" is from 1.57-1.61 /im, near a peak in planetary emission, 
as shown in Figure |T0] The state of the planet at the end of the core-accretion method of formation, from lHubickyj etaT] (20051) . is highly 
uncertain, so these predictions should used with care. See .Marley et al.. i .2007b.) for a discussion. 
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2.4 F 




10° 10' 10^ 10' 10° 10' 10^ 10' 



Time (years) 

Fig. 1 . — Planetary thermal evolution models, updated from'Marley et al.' (2007b). Dotted lines indicate "hot start" planets with an arbitrary initial condition. Solid 
lines indicate planets with an initial model from the Hubickyj et al. (2005) core accretion formation inodel. The model atmosphere grid is Ix solar and includes 
the opacity of refractory cloud species. As in lMarlev et aL t,20Q7h> . times on the x-axis are years since formation, which takes no time (by definition) for hot start 
planets, and '^2. 3-3.0 Myr for core accretion planets. 
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500 1000 1500 2000 2500 

T (K) 

Fig. 2. — Cloud-free P-T profiles at 1400, 1000, 600 K at log g=3.67. Curves in black are for 1 X solar metallicity. Curves in red are for 5x solar metallicity. 
Filled circles indicate the pressure of the mean photosphere, where T=T^ff. Dotted curves show locations of cloud condensation while dashed curves are chemical 
equal-abundance boundaries. Only the 1 X boundaries (black) are labeled. Note that condensation curves shift to higher temperatures as metallicity increases, while 
equal-abundance boundaries shift to lower temperatures. 
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1 3 10 

Wovelength (fim) 



Fig. 3. — Emergent spectra (erg s"' cm"^ Hz"' ) for 5 X solar ([M/H=0.7], colors) and 1 X solar models ([M/H=0.0], black) at, from top to bottom, 1400, 1000, 
700, and 500 K, for log g=3.67. The inset shows the 700 K models on a linear x scale (from 0.65 to 5 fim) and linear y scale (x 10"', relative to the rest of the 
figure). 
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1 2 3 4 5 6 

Wavelength (/^m) 

Fig. 4. — Flux density at 10 pc for 4 Mj objects at an age of ~10 Myr. In red is a hot start evolution model with solar metallicity at 1000 K. In light blue is 
this same model with 5 X solar metallicity, for comparison. In dark blue is a 600 K model that uses the core-accretion initial condition and 5 X solar metallicity. 
Over-plotted in black are 10"'' and 10"' contrast ratios relative to two blackbody stars. The two solid curves are for a Sun-like 5770 K star and the dashed curves 
are for an M2V-like 3600 K star. 
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1 2 3 4 5 6 

Wovelength (/^m) 

Fig. 5. — This is the same as Figure|4] on tlie same scale, but now at age of 80 Myr. The hot start model has cooled to ~600 K, and the core accretion model to 
~500 K. 
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T (K) 



Fig. 6.— Solar metallicity P-T profiles at 1200 and 600 K at log g=5.5, 4.5, 3.5 (thick to thin lines). 
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Fig. 7. — Difference in magnitude in a given band ("metallicity color") as a function of T^ff at tliree surface gravities tliat span tlie range of young Jupiter surface 
gravities sliown in Figure[T] Metallicity color is determined by subtracting the magnitude of the 1 X model from the 5 X model. Implicit is the assumption that both 
planets have the same radius. For instance, at log g=3.67 and reff=1000 K, the 5x solar model is redder in J-K by 0.55 (0.25-(-0.30)). 
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Fig. 8. — Absolute magnitudes vs. time for a 4 A/j planet. In solid lines is the 5 X core-accretion start model while in dashed lines is the 1 X arbitrary hot start 
models. 
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Fig. 9. — Emergent spectra (erg s"' cm"^ Hz"' ) for 5 X solar ([M/H=0.7], colors) and 1 X solar models ([M/H=0.0], black) at, from top to bottom, 1400, 1000, 
700, and 500 K, for log g=3.67. Thin lines are for equilibrium chemisti'y, as shown in Figure|3] Thick lines show models that utilize non-equilibrium chemistry with 
log X;;=4. Infrared filter bandpasses are shown in gray on the top panel. 
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Fig. 10. — A four-panel comparison of flux densities (erg s"' cm"^ Hz"') around //-band. All are referenced to a model with Teff=70() K, log g=3.67, solar 
metallicity, and equilibrium chemistry. The T^n panel shows a ±100 K change in T^ff. The [M/H] panel compai'es the standard model to one that is 5x solar 
metallicity, [M/H=0.7]. The upper right panel (log g) shows the effects of gravity, while the lower right panel shows the effects of non-equilibrium chemistry due to 
vertical mixing. Gray vertical lines guide the eye and illustrate possible locations of narrow band filters. 



